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We extend our recent analysis of the dilepton and photon 
emission rates to the case of finite temperature and baryon 
density, within the context of a density expansion. To leading 
order, the effects of the baryon density are assessed using data 
(photon emission) or constraints from broken chiral symmetry 
(dilepton emission). Next to leading order effects are worked 
out, and their contribution qualitatively assessed. The open- 
ing of the ttN cut causes the photon rate to saturate the em- 
pirical photon yield of WA80, but may not be enough to fully 
account for the excess low mass dileptons seen at CERES. 



1. Introduction 

Recent relativistic heavy-ion collisions at CERN have 
reported an excess of dileptons over a broad range of lep- 
ton invariant mass, below the rho mass |,|. A possible 
excess was also reported in the direct photon spectrum 
A rate departure from p-A collisions may indicate 
some medium modifications in the hadronic phase 0K an 
idea that has spurred considerable attention lately |^-^ . 

In a recent paper (here after referred to as I), we 
have analyzed the dilepton and photon emission rates 
from a baryon free hadronic gas using on-shell chiral re- 
duction formulas [ pT| in the context of a density expan- 
sion. At pion density n,r and temperatures of the order 
of the pion mass m^, the expansion parameter was iden- 
tified with = n^/2m,r/^ J; 0.3 [0, where /^r = 93 
MeV is the pion decay constant. Some enhancement 
was observed in the low mass dilepton and photon rates 
due to the "tails" in the vector and axial-vector correla- 
tors. This enhancement is however insufhcient to account 
for the empirical excess when inserted in a conventional 
hydro-evolution of the fire ball, and could be traced back 
to the detector cuts in transverse momentum . To ac- 
count for the data, an order of magnitude enhancement 
in the present bare rate is needed in the low mass region 
from 2m^ to Am^^. 

At present CERN energies, both the S-Au and Pb-Au 
experiments show evidence of a sizable baryonic density 
in the mid-rapidity region (about two to three times nu- 
clear matter density po). Could it be that such an en- 
hancement is caused by the strong pion-nucleon interac- 
tions? Some recent investigations suggest so in the con- 
text of some effective models of pion-nucleon dynamics 
and at specific kinematics (back-to-back emission) |^,^ . 

In this letter we would like to investigate the role of a 
finite baryon density in a hadronic gas using the approach 



presented in I, with an emphasis on constraints brought 
about by broken chiral symmetry and experiments be- 
yond the threshold region. Specifically, we will use the 
comprehensive framework introduced in in the form 
of on-shell chiral reduction formulas. The outcome is 
a set of general on-shell Ward-identities that could ei- 
ther be saturated by experiment, or approximated near 
threshold by an on-shell loop-expansion |ll|,p^ , and 
above threshold by resonance saturation ||14[] . The re- 
sults conform with the strictures of broken chiral symme- 
try, relativistic crossing and unitarity. In this sense they 
are more than conventional dispersion analysis. Both the 
photon and dilepton rates are subjected to the same anal- 
ysis due to their common dynamical origin, albeit with 
different kinematics. 

In the presence of baryons, we are still able to organize 
our calculation into a density expansion as was done in 
I for the pions alone. This allows us to separate out the 
important effects and ascertain the corrections from the 
higher order terms. In section 2, we outline the charac- 
ter of this density expansion at finite temperature and 
density. In section 3, we discuss the photon rates, us- 
ing data and chiral constraints. The effects of nucleons 
are assessed using a relativistic one-loop chiral expansion 
without the A and a chiral expansion with the A is eval- 
uated in section 4. In section 5, we extend our discussion 
to the dilepton emission rates to determine the enhance- 
ment due to nucleons. We then assess the magnitude of 
the next to leading order corrections in the density ex- 
pansion in section 6. Our conclusions are summarized in 
section 7. 

2. Dilepton and Photon Rates 

In a hadronic gas in thermal equilibrium, the rate R 
of dileptons produced in an unit four volume follows 
from the thermal expectation value of the electromag- 
netic current-current correlation function [T^ ] . For mass- 
less leptons with momenta pi,P2, the rate per unit in- 
variant momentum q = pi + P2 is given by 



dR 



W{q) 



(1) 



where a = e^/Air is the fine structure constant, 

W{q) =J d4a;e-**-"Tr(e-("-^^-^^)/^ J''(a;)J^(0)) , (2) 

eJp is the hadronic part of the electromagnetic current, 
H is the hadronic Hamiltonian, /i the baryon chemical 
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potential, N the baryon number operator, the Gibbs 
energy, T the temperature, and the trace is over a com- 
plete set of hadron states. For leptons with mass to/, the 
right-hand side of (0) is multiplied by 



(1 



2mf 



Similarly, the rate for photons follows from (|^) at = 0. 
Specifically, 



^ d^q 47r2 



(3) 



Both (0) and (|) follow from the same current-current 
correlator (||), albeit for off-shell and on-shell photons 
respectively. For consistency, both emission rates will be 
assessed simultaneously as was stressed in I. 

From the spectral representation and symmetry, the 
rates may also be expressed in terms of the absorptive 
part of the time-ordered function 



W(q) 



1 + e9"/T 



ImW^((7) 



(4) 



W^{q) yd^a;e'«-^Tr(e-("-^'^-^'/^ T*J^(a;)J^(0)) . 

We observe that (^) vanishes at T = even for non- 
zero baryon chemical potential. A cold nuclear state can 
neither emit real nor virtual photons, for otherwise it 
would be unstable. 

For temperatures T < m^r and baryonic densities 
riN ~ 3yOo we may expand the trace in (|^) using pion 
and nucleon statesQ Expanding the trace in terms of 
free pions and nucleons with Feynman boundary condi- 
tions, and summing over disconnected pieces lead to the 
density expansion 

~iW^{q) - (0| Oj(g)|0)eonn. + l '^^ (?) K) conn. 

+ J2 ^^j(9)l^>conn. + ... (5) 

s,I 

with the nucleon fields carrying implicit spin and isospin 
N — {p, s), the pion fields implicit isospin tt = 7r''(fc), 
and 

The phase space factors are 



'^We neglect the anti-nucleon contribution since it is highly 
suppressed for the temperatures and densities we consider. 
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(27r)3 2Ep e(^p-^)/r + 1 
d^k 1 1 
(27r)3 2ujk e'^"/^ - 1 



with the nucleon energy Ep — \J mj^ + p"^ and the pion 
energy = y'ro^H-A?. The matrix elements in (||) cor- 
respond to the forward scattering amplitudes of a real 
(photon rate) or virtual (dilepton rate) photon with on- 
shell nucleons and pions. We consider the nucleons to be 
in thermal equilibrium and ignore the off-equilibrium ef- 
fects caused by the collision dynamics on the mid-rapidity 
nucleons. The chemical potential is fixed by specifying 
the nucleon density. We observe that the first term in (J^) 
when used in conjunction wit h (H ) is reminiscent of the 
resonance gas approximation |lq |. The density expan- 
sion is an expansion in the number of pions and nucleons 
in the final state, which means by detailed balance and 
time-reversal all possible thermal reactions in the initial 
state. 

To leading order in the nucleon and pion densities 
= 4 J dN 2Ep and rij^ = 3 J dir 2wk only the first 
few terms in (||) contribute. A slight refinement is to note 
that the reduction of one pion is associated with a factor 
of I/Ztt, and that of a nucleon is associated with a factor 
of qa/ f-m where qa = 1.26 is the nucleon axial charge. 
The dimensionless expansion parameters should be = 
n7r/2m7r/^ and kn = njs[ g\/2mj\j f^. For T < and 
riN < 3poj we have ^ kjv ^ 0.3. The truncation 
should be reasonable, unless severe infrared divergences 
develop. This is unlikely in the hadronic phase since most 
resonances are massive with finite widths. 

3. Adding Nucleons 

A full treatment of the first two terms of eq. (||) was 
carried out in I. There the electromagnetic current was 
decomposed into an isovector part and an isoscalar 
part B. The BB and BV correlators are expected to be 
small in the pionic states, but in the presence of nucle- 
ons they must be reassessed. The vacuum contribution is 
proportional to Ily, the transverse part of the vector cor- 
relator (0|T*VV|0) which follows from electroproduction 
data. Also including terms to first order in the density 
expansion, we find 

ImW^(g) = -Sg^Im Uviq^) 

+ J_ JdTTW^{q,k)+ ldNWf,{q,p) 

+ 0{kI,k%,k.^kn) . (6) 

The term linear in pion density can be reduced by the 
use of chiral reduction formulas to a form amenable to 
experimental determinations. That analysis showed that, 
to within a few percent, the only important contributions 
were 

Wf(g,fc)~ 12921m Uy{q') 
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-6{k + q)Hm ((fc + + (g - 
+ 8((fc • g)2 - m2q2)ini Uv{q^) 
xRe(Afl(fc + (j) + Afl,(fc-g)) 
with Afl(fc) the retarded pion propagator 
1 



(7) 



A^(fc) = PP 



fc2 



-Z7rsgn(fc")(5(fc2-m^), 



and TIa the transverse part of the axial correlator 
(0|T*jyijA|0) which follows from tau decay data. 
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FIG. 1. The photon emission rate to first order in the 
density expansion for pions only (dashed) and both pions and 
nucleons (solid). A fixed nucleon density of po was used. 

Now taking into account the nucleons gives the third 
term, which is just the spin-averaged forward Compton 
scattering amplitude on the nucleon with virtual pho- 
tons. This is only measured for various values of < 0. 
However, the dilepton and photon rates require > 0. 
Therefore, only the photon rate for this term can be de- 
termined directly from data by use of the optical theorem 



(8) 



with s — {p + qY . Experimentally, the isospin sum of the 
cross section in (||) can just be replaced by cr^ot. [0 '^^^ 
to the weak collective effect of the two nucleons in the 
deuteron. The total cross section of 70? is dominated by 
the A(1232) The result for the instantaneous photon 
rate is shown in fig. ^ for T = 50, 100, 150 MeV and a 
nucleon density of njv = po • From there it is easy to see 
that the nucleon contribution dominates over the pions 
by about an order of magnitude. If the same situation 
were true for the dilepton rate, the enhancement would 
probably be enough to explain the excess dileptons in the 
low-mass region. 



This huge correction to the pion result near threshold 
is due to the opening of the A channel. Approximating 
the cross section in (H) by a narrow A resonance, we 
see the dimensionless parameter of importance is (toa — 
mAr)/2rA ^ 1 which is multiplied by k^t under thermal 
integration. This enhancement from the opening of a new 
threshold does not iterate in the low mass region. This 
ensures the density expansion is under control. In section 
5 we check whether this enhancement in the photon rate 
is enough to surpass the upper bound set by WA8G |^ . 

For off-shell photons, we must resort to chiral con- 
straints to determine the nucleon contribution to the 
dilepton rate. Broken chiral symmetry dictates uniquely 
the form of the strong interaction Lagrangian (at tree 
level) for spin i. Perturbative unitarity follows from 
an on-shell loop-expansion in 1/ Jtt, that enforces current 
conservation and crossing symmetry |]l3| . 

In general, after summing over spin, only four invariant 
functions which depend on s and q^ remain: 

/ d4xe'^--(7Voutb)|r*J^(a;)J,(0)|7Vin(p)) 

s 

= 4 {g,,„A + q^quB + q(f,Pi,)C + Pf,p^D) . (9) 

The absorptive part, which is of relevance here, only 
starts to contribute at one-loop. The four invariant func- 
tions are related by current conservation — leaving only 
A and B as independent. To one-loop in I/Ztt the dia- 
grams that contribute are shown in fig. |^ and require the 
addition of the crossed diagrams as well. The analytic 
expressions are quoted in Appendix A for arbitrary q^. 
For on-shell photons, they fulfill the non-rcnormalization 
condition A{rrij^,Q) — Z with Z — 0(1) for the neutron 
(proton). This is fully given by the tree level nucleon 
contribution, leaving ^'°°P*'(m^, 0) = which has been 
checked to hold analytically for the loop functions. 
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FIG. 2. One-loop diagrams for Compton scattering. The 
graphs with a star also have a mirror image diagram which 
must be taken into account and all graphs except for those in 
the last line require the addition of a crossed diagram. The 
lettered graphs are discussed in Appendix A. 



The results for the photon rate are shown in fig. g 
for a representative temperature T = 150 MeV (results 
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are similar for other temperatures). Comparison is made 
with the nucleon result from fig. The one-loop ampli- 
tude reproduces the shape attained from the data quite 
well — even for larger photon energies — with the ex- 
ception of the small bump around the A peak. That 
missing feature, however, is kinematically around the re- 
gion of interest for the enhancement of the dilepton data 
(300 MeV above threshold). Therefore to make a full 
statement about the nucleon contribution for finite , 
we will include the A in the analysis. 




FIG. 3. The photon emission rate as in fig. 1 using 
data and the optical theorem (solid), using the one- loop 
results (dashed), and using a contribution from the A 
(dashed-dotted). A fixed nucleon density of po was used. The 
inset shows data for cr^ot in millibarns as a function of lab 
momentum in GeV from [17] and the good agreement for the 
imaginary parts of the one-loop and delta results discussed in 
the text. 



4. The Role of the A 

Chiral symmetry does not constrain the form of the 
strong interaction Lagrangian at tree level for spin | in 
the relativistic approximation. We therefore treat the A 
as a resonance in the continuum, and use time reversal 
invariance and data to constrain the form of the tran- 
sition matrix element. We decompose the TV- A vector 
transition matrix element in general as^ 



■^We could also include terms with 7^u'', but these are sup- 
pressed near the mass-shell, where we are interested. 



{N{p)\Vl\^{k)) =u{p) 



{R{q^)+4R(q^)) {q^q,. - g^.q'') 



iS{q )cr^qxq^ 



(10) 



with q = k — p. Current conservation dictates the spe- 
cific tensor form of ([lo|). Note that the isoscalar part of 
the electromagnetic current does not contribute to the 
nucleon-A transition matrix element. 

The modulus square and sum over all spins and 
isospins of ( |l0| ) will be denoted by M.{k'^,q^) and is 
quoted in full in Appendix B. It is directly related to 
the A contribution of the forward scattering amplitude 



W^(g,p)=Im- 



M{s,q^) + {s 



u) 



with u = [p — qY . If we were not to take the imaginary 
part, the full tree level result must also include the nu- 
cleon Born terms. At = this is equal to the isospin 
sum of 8A(s,0) and the same non-renormalizability con- 
dition of the last section requires A'^(m^, 0) — 0. This is 
satisfied trivially by the decomposition used in eq. (p^), 
since the charge was fixed to zero by current conserva- 
tion. 

If we do not contract the vector indices of as in 
eq. (H), we can also determine B{s,q^). The polarizabil- 
ities are defined as 



a 



+ (i = -2amNA"{m%,0) 



a 

TUN 



P = B{ml,Q) 



8a 



n2 



TO/ 



'■N 
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over-constraining the value of Q. Although experiment 
has shown the electric polarizability dominates over the 
magnetic one a > /3, the A is a magnetic dipole effect 
and tends to overestimate (i in other calculations p9| . 
Therefore we concentrate on fixing a reasonable value 
for a + /3 = 10 X 10"^ fm'^. This indeed gives a ~ 
14 X 10""' fm'^ about three times the experimental value. 
This corresponds to Q(0) = 2.75/to7v. Translating into 
our notation, it is comparable to the value Q ~ l.hjmj^ 
obtained by a different method [2C| ]. 
The decay width of the A is given by 



rA^Ar7 = aiTiN- 



8m\ 



7W(TOi,0) ~ 0.72MeV, 



1.2/m^. 



This compares 



and can be used to fix S{0) 
well with 5" ~ 1.0/m^ from 

The form factors R, R do not contribute at q^ = 0. For 
finite g^, the A contribution is mostly around its mass 
shell and so we can use the Dirac equation to rewrite 
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(mA + mN)R in ( p^ and absorb R into the defi- 
nition of R. Then using the vertex ( |lo|) in pion electro- 
production, R can be shown to be proportional to the 
longitudinal part of the cross section. This is on the or- 
der of 10% of the total result even for moderate and 
therefore we will take R = R ^ for the rest of this 
paper. 
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FIG. 4. The dielectron rate for pions alone (dotted), pions 
and A (solid), and pions and one- loop (dashed). The contri- 
bution from pions. A, and one-loop together is represented by 
the thick solid line. A fixed nucleon density of po was used. 

We plot the photon rate due to the A along with the 
one-loop result in fig. ||. From there we see that the A 
rate does quite well around the peak and then dies away 
for larger energies. In fact, adding the one-loop and A 
rates, we come very close to reproducing the nucleon re- 
sult obtained from the data. This even includes the sec- 
ondary peak, seen in the inset of fig. ^, which is fit mod- 
erately well by the loop contribution. Therefore we take 
both contributions into account as a way to parameterize 
the photon spectral function with on-shell chiral symme- 
try, relativistic crossing, and unitarity constraining the 
function near threshold, and experiment constraining it 
above. We now can move to finite to ascertain the 
dilepton rate. 



5. Dilepton Rates with Nucleons 

We will neglect the momentum dependence of the 
form factors in eq. ( |lO| ) since they are expected to be of 
monopole form {1 — q^iM^)^^ with A/+ = i(mjv + mA) 
and hardly change for the q^ we will examine (less 
than 1 GeV^). Then taking the results of the previous 
section to > gives the curves in fig. ^. Also shown 
is the one-loop result and the pion result from I for com- 
parison. 

Qualitatively, the rates are similar to the photon case. 
The enhanced bump of the A is smeared out and en- 
hances the rate by as much as a factor of ten (for 
T = 100 MeV) before gradually dying out again near 
the p peak. The tail of the one-loop goes above the A 
result around M = 200 MeV. This may be easily under- 
stood by noting that the absorptive part involves the irN 
cut, with pions typically carrying a momentum of order 
rriTr and the nucleon a momentum of order pp ^ m^. 
The one-loop is also overtaken by the pion result at 
the p peak. A closer examination of the one-loop real 
part shows that it is about half of the tree result for 
q^ ^ (SOOMeV)^, above which the l/f-K expansion is no 
longer reliable. Fortunately, the energetic dilepton pairs 
with g" > 300 MeV are thermally suppressed by almost 
an order of magnitude at the highest temperature, typi- 
cally e-300/150 ^ 1/7. 
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FIG. 5. The dielectron rate including the pions. A, and 
one- loop for po (solid), ipo (lower dashed), and 3po (upper 
dashed) . 

Combining the one-loop and A results gives the thick 
solid line in fig. IJ. The dominant effect in our case 
comes from the continuum and not the A resonance. At 
M = 400 MeV, the inclusion of nucleons enhances the 
rate by a factor of three. Others who have taken nu- 
cleons into account through various methods |§,||J|] find 
enhancements in the rate similar to our result. However, 
the details of the in-medium effects are not needed in this 
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analysis and so should be considered more general. Also 
in our calculations there is no shift of the dilepton pair 
production p peak. 

In fig. ^ we show the dependence of the T = 150 MeV 
dielectron rate on the nucleon density for = ^pa, po, 
and 3po- At T — 150 MeV and ^pq, the result is dom- 
inated by the pions as the nucleon contribution is ex- 
tremely small. This shows the effect is truly from the 
nucleon density. Even so, only at 3po does the rate flat- 
ten out at the value of the p peak rather than increasing 
towards it. 

In order to fully understand the role of the cuts, we fol- 
low ^ and evolve the instantaneous rate over space-time 
using recent transport equation results |^ and experi- 
mental cuts [|l|. Assuming a homogeneous expansion, 
the volume and temperature only depend on time and 
can be parameterized as 

V{t) ^Vo(^l + T{t) = (T, - T^) e-'l- + 

with to = 10(10.8) fm/c and r = 8(10) fm/c for S-Au 
(Pb-Au) collisions. The freeze-out time is tf.o. = 10(20) 
fm/c. We absorb Vq into an overall normalization con- 
stant that includes the charged particle distribution as 
well. RQMD predicts the initial baryon density ub — 
2.5po for S-Au and 4po for Pb-Au ||2^ which translates 
into nN — 0.7po(1.0po)- Only the pion and nucleon den- 
sities are retained since we are expanding the emission 
rates (|) in terms of stable final states, summing over all 
possible initial states. This point is particularly trans- 
parent in the p region. Assuming chemical equilibrium 
gives a constant nucleon chemical potential. The rate 
then can be written as 



{dN/dri dM) 
{dN,h/dv) 



= NqM 



dtVit) /^A(,«,,^)^ 



with the normalization constant A^o = 6.76(1.33) x 10^^ 
fixed by the transport results in The acceptance 
function A ensures the detector cuts at CERES (p± > 
200(175) MeV, 2.1 < r; < 2.65, and O^e > 35 mrad) 
are taken into account. The finite mass resolution is also 
taken into account by folding the spectrum with a Gaus- 
sian averaging function given by the CERES collabora- 
tion |l|l. 

Finally, we plot the results for S-Au and Pb-Au colli- 
sions in fig. ^. This requires the addition of the Dalitz 
and omega decays which were obtained from the trans- 
port model Q . The data is presented by adding the sta- 
tistical and systematic errors linearly. The Pb-Au data 
is still preliminary. For hn — 0, the pion result alone has 
already been analyzed in two hydro-dynamical models 
[p^,|6| . Comparison of this alone can be used to test the 
validity of the evolution in space-time described above. 
Our result is shown as the dashed line in fig. ^. The over- 
all shape due to the cut looks reasonable. The values for 
the rate in p^] are an order of magnitude smaller, but 



this can be traced back to a shorter time in the hadronic 
phase. The agreement with [|| are fair. We can think of 
the results of this space-time evolution as an upper bound 
on the results from a true hydro-dynamical model. 

Adding the nucleon contribution gives the solid line 
in fig. ^. The effect of the cuts is dramatic, resulting 
in a very small enhancement. Only if we take the ex- 
treme case of the baryon density totally saturated by 
nucleons do we start to reach the lower error bars of the 
data in the M = 200 — 400 MeV regime as shown by 
the dashed-dotted line. The large effect seen in fig. ^ is 
not present because the temperature dies away quickly, 
thereby decreasing the nucleon density and rate dramat- 
ically as seen in fig. ||. 
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FIG. 6. Our dielectron rate including the A and one-loop 
contributions evolved in space-time as in [8] for S-Au and 
Pb-Au collisions. In the upper graph, njv = 0,0. 7po, and 
2.5po are plotted as the dashed, solid, and dashed-dotted 
lines respectively. In the lower graph, the lines are for 
njv = 0, l.Opo, and 4po. The data are from [1]. The sys- 
tematic errors are added linearly to the statistical error bars 
to give the cross line. The Pb-Au data is preliminary. 

Comparing to our rate is slightly smaller in the 
region of interest. This can be traced back to the bare 
rates. For M > 400 MeV, their rate for tin — 0.5po is 
larger than the rate for hn = po [ p3| which differs from 
our result in which the nucleon contribution dies away 
continuously leaving only the pion contribution (fig. ||). 
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This is related to the depletion of the p peak in [g| , which 
does not occur in our case. 

We can also evolve the photon rates and compare with 
the upper bounds set by WA80 for S-Au The pro- 
cedure is similar to that for the dileptons above. The 
acceptance cut is 2.1 < rj < 2.9 and no integration over 
momentum is required. In this case the factor TVq is only 
the volume Vq which we normalize by comparison with 
[ p2[ for riAT = 0. This fixes Vb — 300 fm^ which corre- 
sponds to an initial radius of 4 fm and agrees with central 
collision estimates. The result for the pion (dashed line) 
and nucleons for the extreme case of n^r = 2.5po (solid 
line) are shown in fig. ^ The inclusion of nucleons put 
the rate right on the edge of the upper limit for the data. 
Since we have analyzed both the dilepton and photon 
rates simultaneously, this implies that more enhancement 
of the dilepton rate would overshoot the photon data, a 
particularly important point in our analysis. 




FIG. 7. Our photon rate with (solid) and without 
(dashed) the nucleon contribution discussed in this paper. 
The data are upper bounds from [3]. 



6. Higher Order Terms 

We have now assessed all the terms of eq. ^ linear 
in the density. The result is that they cannot explain 
the CERES data. In general, higher order terms will be 
suppressed by more powers of the expansion parameter. 
This implies the inclusion of higher order effects will not 
be dramatic. However, this rule of thumb can be upset 
if new thresholds open up in the low mass region. Typ- 
ically, the second order corrections to eq. ^ are of the 
form 

<5ImW^(g)~^ jd^^d^2^^,{qMM) 



1 

T2 

J 7T 



2 , d7rdNW^^iq,k,p) 



i JdNidN2W^^{q,pi,P2) 



(11) 



with W^^ referring to the forward scattering for real 
or virual photons ^*a(3 7* a/? as used above. Note 
that matrix elements involving nucleons carry an implicit 
factor of 1//^ since the contribution only starts at one- 
loop. Experimentally, we can estimate the strength of 
the third term in eq. ( [ll|) at = as follows. We noted 
in connection with eq. (8) in the last section that the 
cross section for Compton scattering on a deuteron is 
identical to the sum of Compton scattering on a proton 
and neutron individually. This seems to indicate that 
the collective effect of two nucleons is small and we will 
ignore it for the rest of this paper. 

The A may be thought of as a bound state of a pion 
and a nucleon, and so the second term in eq. (^Tj) could 
be important. We also must retain the isoscalar contri- 
bution to the electromagnetic current when nucleons 
are involved. We can still reduce out the pious as given 
in ||ll[] if we assume the pious do not interact with the 
isoscalar part. This is true in the soft pion limit and so 
should be a good approximation [ p^ . The overall form 
of is very similar to W^, in which one pion was 

also reduced out, with the vacuum matrix elements being 
replaced by averages in a nucleon state. The result is 



4^ ^iq-x 



W^^(g,fc,p) = -4Im*^ Jd 

S,I 

x(A^out(p)|T*J^(a;)J^(0) - B^(a;)B^(0)|Mn(p)) 

_|_ gaSg^aS/ 
a,{s,/} 

xlmf 5^-^ - (fc + g)(^(2fc + qY^^R{k + q) 



+ {k + qT{k + qY{2k + qfA\{k + q) 



X j d'x e^('=+'')-z(7Vout(p)ini(a;)jL(0)|Mn(p)) 
+ (9 ^ -q) 

s,I 

xRe{ARik + q) + AR{k-q)) 
x{NoM\T*J^{x)3A0) - B^(a;)B,(0)|M„(p)) 



J2 Im (iVout ip)\T*3^ (x) (y)<7(0) I A^in ip)) 



S,I 



-k^'kl^ Jd'^xd^yd^ze"'-^''-y'>e-'''-'' 

X Im z ( Aout ip)\T*3^ {x)3^ (y)jl„ (z)jl^ (0) | A^ {p)) 
+ kf'lm (<5^ - (2k + q)^{k + qTAnik + q)) 
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X ^ ^6"3^(iV„,,(p)|T*j^Jx)jl^(y)J^(0)|Mnb)) 

a,{sj} 

+ {q^ -q) + {k^ -k) + iq,k^ -q, -k) 

with Ua denoting the transverse part of (A^|T*j^j^|A^) 
summed over spin and isospin. We use the same analysis 
as the last section for the JJ term. The contribution 
here is negative, but weighted by an extra factor of n. 
It is about 10% of the leading contribution. Since there 
is no prominent resonance with the quantum numbers of 
the isoscalar term BB, the imaginary part should also be 
small and will be neglected. 

The second term contains the axial current j^, which 
is a pionic contribution, and so this term may feed into 
the A. Using the decomposition (k — p' — p) 



)kf^ki. 



Gik^h^k, 



il{k^)a^kxk^ 



an analysis of ttN scattering in [gj] showed that F[Q) = 
1.38 and G(0) = 0.4/mAr. Using the fact that H and / 
are suppressed for large iVc, we can justify keeping only 
the F term to find 



9 " {p + Ky-m\+ iniATA 
{p + K)^{p + K), 



[K -K) 



where only the A channel has been retained. This term 
and its crossed counterpart gives a positive contribution 
to the rate. Numerical results show this term is apprecia- 
ble right near threshold but quickly shrinks to 10% of the 
terms linear in k, and practically vanishes even before the 
A peak is reached. The damping is simply coming from 
the exponential factors in phase space and so the only 
place where this term is significant is within the Dalitz 
tail. This behavior is expected to occur for all terms in 
and so we may neglect the terms of order kjvKtt. 
For completeness, we quote the result for the first term 
in (|l|) in Appendix C. Qualitatively, we note that the 
full TTTT scattering amplitude as well as terms of photon- 
pion scattering similar to Wtt appear with an additional 
suppression factor of k. In the soft pion limit most of the 
correlation functions in the pionic state are amenable to 
correlations in the vacuum, some of which were assessed 
in I and found to be small. Hence, we expect qualitatively 
that W^^ is suppressed by the factor as expected and 
hence are small in comparison with the lower order terms. 



7. Conclusions 

We have extended the analysis in I to account for pos- 
sible nucleon densities in assessing the photon and emis- 
sion rates from a hot hadronic gas. We have used the very 
general framework put forward in jl^Jl^,|l^ to allow for 
a comprehensive analysis of the emission rates in pow- 
ers of the matter densities, and enforce the known con- 
straints of broken chiral symmetry on-shell, current con- 
servation, and relativistic unitarity. We believe that most 
approaches to this problem as well others involving issues 
of chiral symmetry in matter, have to embody the chiral 
Ward- identities established on-shell in [|ll|, whether at 
threshold or above. 

Our purpose of this paper was to use these arguments 
in combination with a density expansion to account for 
a systematic analysis of the dilepton and photon rates 
from a hot hadronic gas at finite nucleon density. For 
the photon rate, we have used the constraints from bro- 
ken chiral symmetry and data to leading order in the 
densities. A large enhancement of photons in the low en- 
ergy region was found. Using an on-shell expansion to 
one- loop and a general decomposition of the A, we were 
able to reproduce the photon rates at = 0. 

We have used the same amplitude to estimate the 
dilepton production rate at q^ > 0. In the latter, we 
have found an enhancement from the ■nN cut in the low 
mass region for a nucleon density on the order of nu- 
clear matter density on top of the already large enhance- 
ment from thermal reactions alone. This enhancement is 
caused by the opening of a threshold in the TriV channel, 
and does not iterate coherently in the next order correc- 
tions, which we have found to be qualitatively small and 
mostly confined to the Dalitz region. 

For S-Au and Pb-Au reactions, we have evolved the 
dilepton rate in space-time and taken the experimental 
cuts. The result is not enough to explain the data, even 
with some liberal choices of the parameters. The nu- 
cleon densities do trigger large dilepton yields but die 
out quickly after short reaction times. In addition, we 
have evolved the photon rates and found the results on 
the edge of the upper bounds of the data in the extreme 
case. This means that any further enhancement would 
overshoot the empirical measurements. Since we have a 
controlled expansion scheme, these results lead us to sus- 
pect the nucleons may not explain the data. Any mecha- 
nism which could explain the data would need to have a 
very specific form in order not to upset the photon yield. 
However, this space-time evolution is still only a qualita- 
tive assessment. Quantitative results from a full hydro- 
dynamical model which takes into account the quark- 
gluon phase, flow, and hadronic multiplicities is needed 
for a final analysis. 

Our results imply a different interpretation for the role 
of the nucleons compared with a number of recent esti- 
mates [pP,p5t. Our analysis, however, is different in a 
number of ways. We have evaluated the rates at finite 
temperature, since there is no emission in cold matter. 
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Also our results go smoothly into those of I at zero nu- 
cleon density, thereby enforcing the known constraints 
of (on-shell) broken chiral symmetry, current conserva- 
tion, relativistic crossing and unitarity for all dilepton 
kinematics. It is important to stress that our density 
expansion is an expansion in the number of pions and 
nucleons in the final state, which means by detailed bal- 
ance and time-reversal all possible thermal reactions in 
the initial state. Finally, our qualitative assessment of 
the next to leading order shows that the expansion is re- 
liable in the low mass region, with no apparent need for 
coherent resummations. 

The empirical consequence of these differences is that 
the dilepton strength in the p region from a hot hadronic 
gas remains about constant and dominated by thermal 
TTTT pairs. The enhancement below this region in invari- 
ant mass does not come at the expense of this peak. A 
consequence of this is the enhancement is very sensitive 
to the temperature and nucleon density and dies away 
quickly as the fireball evolves. The resolution in the p 
region for the CERES data is not enough to sort out 
the dilepton emission from 'cold' w's. An improvement 
on this resolution maybe an important step in assessing 
the effects of interaction as reported here in comparison 
to those discussed in [^,||j|j2^. Another improvement, 
as stressed in I, is to have more empirical correlations 
between the photon and dilepton yields. 
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Appendix A 

The on-shell one-loop result can be written in terms 
of Feynman parameter integrals. Since we are only con- 
sidering the imaginary parts, the (real) subtraction con- 
stants will be ignored. 



flA{s,q^)^Z{gl-l) [l- 



'TV 



'■N 



(s - m\Y 



G-'[l-^Z]{s-m%)r-^{s) 



2ZG\s - m^)r^^(s) + 2ZG'q'Tf'' (s) 



'2J2-r--n:N I 



[Z - 2)2miG'gr +{Z~ l)AmiG'ni 
AmlG^g^^ + \{Z - 2)G\^ [mlg^^ + 2rf/(q2)] 
crossed 



flB{s, q^) = -4 (mjvCT^Ar - {gAruN - Gf) 



ZG-^ 



'AT 



2rj^(s)-K2rf"(s)-r'"'^(s) 



N I 



+ 4rf "(s) -I- 2r^^'^(s) - AT^is) - 4rf "(s) 

2{z-2)G' [rf5^(9^)-rr(<z^)] 

2ZG''Tlg{s) ~ 2G''Tl^{s) ~ G^ [r^'^(s) - 2rf/(s)] 

2ZG^ re(.)-rf/(.)] 

2{Z - 2)mlG^ [gi^ - ^3^^] + 4(Z - l)mlG^il^ 
1 



AmlG"^ 
crossed 



4^ 



c = 



q^ 



'AT 



{A + q^B) D = 



2q^ 



C 



m 



N 



with r456 = r4 + 2r5 + Tq and Fia = Fi -[^ Fs. Also 
gA is the nucleon axial charge, ct^at is the pion-nucleon 
sigma term, G = fTrgnNN, and G = 2G — gAiriN- The 
values used were gA = 1.26, niN — 940 MeV, = 140 
MeV, and ct^tv — 45 MeV to ensure the nucleon is on- 
shell. Taking q^ — and a^^N reduces the loop 
result to that given in |19|. Conformity with the Ward- 
identities |ll| requires an on-shell expansion |13| as taken 
into account here. 

Each of the J, F, g and Si's are Feynman parameteriza- 
tion integrals quoted in full in jist . The imaginary parts 
have cuts for s > sq = {itin + Wtt)^ and q^ > 4m^. The 
cuts for q^ > 4m^ do not impact this paper and are not 
quoted, thereby making the contributions from and 
F^^ zero. The labeled diagrams in fig. |^-h are domi- 
nated by J'^^, V^^, F^'^, g''^, n, g^-^, F'^'^, and J""-^ 
respectively (fully given by these if the nN coupling were 
a constant rather than proportional to a derivative in the 
Lagrangian) . 

We only need to quote a subset of the integrals and 
obtain the others by use of algebraic manipulations as 
outlined in |19]. Using X{x, y, z) — + + — 2xy — 
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2xz - 2yz, K = X{s,m%,ml), Xq = \{s,m%,q'^), and 
Att^ = A((7^, TO^, m^), the imaginary parts for > are 

'A,, 



Im J'^^(g2) 

Im .r^{s) - 
Im r''''{q^) 

Im r'^^(s) = 



Imr^'^(s) 



167rq2 



167rs 
1 



0{q'-^rnl) 
[s - So) 



IGttA 



-4ml) 



7rq 

1 



167r\/X 



:ln 



1 



167r^/A^ 



:ln 



yA^rAg 



(s - So) 
{s - So) 



+ 



1 



167r-v/A^ 



:ln 



Im g'''^ = 

Im = 



All 



^7rg Attq Ag 

S 



_ ) 



wXtt S 
477 - X^Xq 



/A 



7rq 



Sirq^ s — m\ 



X < In 



167r-yA^(s — — q^ 



0{s - so) 

0{s-so) 

0{q^ - 4ml) 
0{s~so) 



+ a/ Att Ag 
~ \/ Att Aq 



In 



+ (s ^ u) 



1 



X < In 



/•rrg + \/ Xj^Xq 
firq — \J X-^Xq 

e{q^ - 4ml) 



^TTQ '^Trg ^q 



In 



mW - \/XqX^q 



with 



a-Kq = {s- m% + + m% - m^) - 2sg^ 

&7rg = (s + - + m% ~ m,l) - 2s{2m% - ml) 

Cirq = q^{s- m% - q^ + 2ml) 

d^q = q (s- mjy - q - m^) 



(s - + q^){s - m% + ml) - 2sq^ 



fyrq = 2s{s - m%) - (s - m% + q^){s - m% + ml). 

These imaginary parts reduce to those in for q^ = 
0. Increasing q^ leads to Xq becoming negative and the 
analyticaUy continued forms of the above expressions, for 
example 



Im r^'^ (s) 



must be used. 



1 



arctan 



\/ XyrXq 



0{s - So), 



Appendix B 

We use the following conventions for the A's sum over 
spins and propagator 



t''t^ \ ^ + m\ 



with 



3 J 2m^'-^'^^^^ 

rVA ij^ + mA) 

3 ) p^-ml + im^T^'^^''^' 



1 1 2 



3mi 



Defining (v, Vp, Vq) — (p ■ q,k ■ p, k ■ q) /mpf and fc^ = s, we 
find 



J2\{Nip)\Vl\Aik))f^Mis,q') 



9mA 



2mNi^qMi{s,q'^) + 2mNvlq^M2{s,q^) 



+ q^Msis, q^) + q^M^is, q^) + q^^M^is, q^) 

Mi{s,q^)^ + Q^ + 2r^VqvS^ 

\ Vq mN J 

- 2r{2v - rVq)QS 
M2{s,q')^(--2r--- 



\] QR 

Vq Vq m\ ) 



+ (vp — m^)R^ - r{rvq + v)RS 

Imi^ 

+ rirvq — v)RR + r^VqVld 



+ - 1 

\Vq\ m 

\ 2 Vq 2toa / 



^ ) — (t'p + "^A) ) RS 



M3is,q') = 



m' 



■ - 2 ) - 3mA ) 



( V ( s \ At 

mi^Vq — — 5 2 + 5r 

\Vq \m^ / m/s 



mA 

f V f s \ 2s 
' mAVq — 9 - 



■v„ 



I QR 



-]-6r+—r)QS 
Vq \ m^ / m/. I 



, f s 2s \ — 

MAs, q~) — Avp — 3mA H n-Vp QR 

\mA m^^ ) 



' (m^Vq — mAi^) (5 ^ I RR 

V "1a/ 

■ ( 3mA - 2fp — ] 

\ mA/ 



'A- 

RS 
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+ v„v brriM r ^(mA + Vpjr K 



'im/^v H 2~(2'7iAri'q — m/\v) ) RS 



M5{s,q^) = - (3 - ^ ) {mA + '^p)lf 

with r = m^/m/y^ and the sum over spins including the 
photon, nucleon, and the A. 

Appendix C 

We quote here for completeness the full on-shell Ward- 
identity for W^j^. It contains the tttt forward scattering 
amplitude and the pion-spin averaged ttj scattering 
amplitude iT-^-y both quoted in 

-^Wf,(q,fcl,fc2) = 

= -(2fci + q)' e^'^e'^'^Im T^'r'^' iih + q), h) 

+ iq^ -q) 

+ 4 {9f.u ~ (2fci + g)^fci.A;j(fci + q)) Im r/^"(q, fca) 
+ {q -q) + ih -ki) + {q, ki -q, -fci) 
+ 1 fc^(2fci + qr J2 ^""^Im (fci, 9) 

+ {q -q) + (^1 -ki) + {q, ki -q, -ki) 



'4y giq-{x~y) 



x{7rUk2)\T*Yl{x)Y^''^{y)am7rtXk2)) 

-72 E ^'^''^e^'^Sfci + g)'^fc?Afi(fci + q)lm B^f^ik,, k2) 

-72 (5^"^ - (^'1 + 9)"(2fci + 9)''A;?Xfci + 9)) 

Jtz 

^Y.e'^^^e'^^nmBliik^ + qM) 

a 

- -^fcf fcf A«(A;i) e^'^e'^'^Im Bl$ik,,k,) 

+ -ifc?A:f (2fci + g)2Afl(fci + g)Afl(fci) 
x5]e'^=^V3^Im6:|(fci,fc2) 

a 

+ (q ^ -q) + (fci ^ -fci) + {q, ki -q, -fci) 

- + + qfC^ki + qf^liki + q) 



-725'^''^? E ^"''I™ + 9, fci) + (fci -fci) 

+Im -^(2fci + g)^(fci + g)"fcf AK(fci + q) 
X^e'^3e^C-Jq,fci+9,fcl) 

a 

-ifc?fcf / d'^xd^yd*z e''=i-(f-^)e^«-^ 
X Im z (^^,, (fc2 ) |T*j^„ {x)rAf3 (2/) (z) V^^'' (0) ki';, (fc2 )) 
with 

^-(fci,9)^^ jd^xd^ye^^^-^e^'^y 
b '' 

xlm (^Lt(fc2)|T* (jlJx)V3(y)) 7r!;,(0)|7rfjfc2)) 
= -^(5'^'=(2fci f d^xe'''^-^ 

Jtt J 

xlm {iTl^,{k2)\T*Yl{x)ami,l{k2)) 

+ Im ^e'^3g(2fc, + 9)^fcf Afl(fci)S:;^(fci, fc2) 

+ ^£^'^'5^ - (2fci + 9)M^f )Im e:«0(fci, fc2) 
-^{ki+qf ( d'^xd'^y e^i ^'e'^' y 

Jtt J 

x\mi{^l^,{k2)\T*Yl{x)rAjy)lApm<{k2)) 

JTT J 

xIm*(7rL(fc2)|r*V^(x)V^(0)|7rf„(fc2)) 
^a?,(fci,M = 

= iYld'x e^'^-^(7rL(fc2)|r*j2iJx)j^^(0)ki'L(fc2)) 

cr-(9.fci+9'fci) = 

= E jd^xd^ ^^(k,+q)■x^-^k,■y 
b 

x(7rL(fc2)|T*V3(0)j5i,(a:)j^Jy)|7rf„(fc2)) 
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